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UNDERV/ATER  EXPLOSION  BUBBLES  III:   THE  EFFECT  OF 
THE  SURFACE  AND  THE  BOTTOM  ON  THE  SHAPE  AND  MOTION  OF  THE 
BUBBLE. 

by  Ignace  I.  Kolodner 

I .   Introduction . 

In  the  absence  of  p;ravlty  and  in  an  unbounded  liquid,  an 
initially  spherical  underi-:ater  explosion  bubble  would  remain  spheri- 
cal, at  all  times.   It  would  perform  radial  oscillations,  v;hlle  the 
center  of  the  bubble  v/ould  remain  at  rest,  and  if  the  liquid  were 
incompressible  these  oscillations  would  be  periodic  and  undamped. 
In  our  first  report,  "Underwater  Explosion  Bubbles  I",  [1],  we 
considered  a  compressible  liquid  and  showed  that  the  oscillations 
would  then  be  damped.   In  our  second  report,  "Underwater  Explosion 
Bubbles  II",  [2],  vie   considered  gravity  and  showed  how  it  x-j-ould 
effect  the  shape  of  the  bubble  and  also  cause  It  to  rise.   In  the 
present  report  vie   considei''  the  effects  of  the  water  surface  and 
bottom,  in  addition  to  gravity,  and  show  how  they  effect  the  shape 
and  motion  of  the  bubble. 

Of  course  all  of  these  effects  have  been  considered  to  some 
extent  by  other  atithors.   However  v:e  have  attempted  in  each  case  to 
present  a  more  complete  and  more  systematic  analysis  than  has  pre- 
viously been  given.   Thus  in  our  first  report  we  showed  that  radiation 
of  energy  from  the  bubble  can  account  for  most  of  observed  energy 
losses.   In  our  second  report  we  showed  that  th';  observed  changes  of 
shape  could  be  accoiinted  for,  at  least  qualitatively,  by  considering 
gravity,  and  that  these  changes  also  diminished  the  rate  of  rise 
and  modified  the  bubble  period  in  the  observed  manner. 

In  the  present  report  in  considering  the  effects  of  the  top 
surface,  the  bottom,  and  gravity,  vre   do  not  constrain  the  bubble  to 
be  spherical  at  all  times,  as  was  done  in  the  previous  work  on  these 
effects  by  Shiffm.an  and  Friedman  [3],  axid  Friedman  [ij.].  We  find  that 
certain  effects,  vrhich  we  call  effects  of  lowest  order,  agree  with 
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those  calculated  by  the  above  cited  authors,  but  higher  order  effects 
also  occur  which  were  not  predicted  by  them.   In  particular  the  change 
of  shape  and  its  consequences  were  obviously  excluded  from  their  work. 
Insofar  as  our  results  coincide  with  the  previous  results,  we  believe 
that  they  constitute  a  proof  of  the  correctness  of  the  previous  re- 
sults, since  our  results  are  based  upon  a  more  accurate  theory.   Of 
course  in  addition  we  obtain  new  restilts  involving  the  change  of 
shape. 

The  present  results  are  based  upon  the  linearized  boundary 
condition  at  the  top  surface  of  the  liquid,  rather  than  upon  the 
exact  nonlinear  boundary  condition.   This  linearization  is  employed 
because  it  simplifies  the  calculation,  and  because  it  does  not  affect 
the  bubble  shape  or  motion  to  the  order  of  approximation  considered 
in  this  report.   Nevertheless  we  are  at  present  considering  this 
problem,  using  the  exact  boundary  condition.   The  results  of  this 
analysis  will  be  presented  in  a  future  report.   They  should  verify 
the  present  results  to  the  order  considered  here,  and  differ  in 
higher  order. 

A  discussion  of  results  will  be  found  in  Section  VII.   This 
discussion  is  qualitative  to  a  large  extent.   Nvunerical  resialts  are 
published  in  a  companion  report,  "Underwater  Explosion  Bubbles  IV  - 
Stimmary  and  Kuraerical  Results",  nO-i-NYU  Report  No.  233. 
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II .   Formulation. 

We  assume  that  an  incompressible,  Invlscld  fluid  bovinded  below 
by  the  rigid  surface   z  =  -  H  and  above  by  the  free  surface 
z  =  z   contains  a  gas  bubble  within  it.   See  Figure  1.   The  velocity 
u(x,y,z,t)   in  the  fluid  is  assumed  derivable  from  a  velocitj'-  poten- 
tial  (j)(x,y,z,t)  v;hich  then  satis- 
fies the  Laplace' s  equation 


z  > 

s 

bubble      /^ 

^ 

^, 

D 

k^- 

X 
\ 

i 

,  J, 

(2.1)  /^  \>  =   0    , 


u 


V  ^>  . 


The  pressure  p(x,y,z,t)   in  the 
v:ater  is  then  given  by  the  Bernouilli 
equation. 


(2.2)   P  =  Pq  -  p[^>t  •*■  |(V  ^)^1-   Pgz 


Here,   P  =  p  +  pgz   ,  v;here  p   is  the  pressure  of  air  above  the 
'    o    o   "^   o  o 

free  surface,   p   is  the  density  of  the  fluid  (water),  and  g  is  the 
acceleration  of  gravity.   Thus,   P   is  the  static  pressure  at  the 
level   z  =  0  . 

The  equation  of  the  bubble' s  surface  is  assumed  to  be  of  the 
form 


{2.3) 


P(x,y,z,t)  =  0  . 


On  this  surface,  the  kinematic  and  the  dynamic  conditions  must  be 
satisfied,  namely 


(2.I4.) 
(2.5) 


Vf  •  V^  +  F^  =  0 


=  P,  -  pih   +  ^(V^)^]  -  P5Z  =  KV'V 


(on  F  =  0 


Thc  last  condition  follows  from  the  assimption  that  within 
tho  bubble  the  prcosuro  is  unif om  and  is  rolatod  to  tho  bubble 
volume  V  by  the  adlabatic  relation,  pressure  =  KV  "^ ,   v;ith  K  and 
Y  constants. 

At  the  bottom  we  have  no  normal  flow,  so 

(2.6)  <|5,(x,y,  -H,t)  =  0  . 

In  a  precise  formulation,  the  top  surface  must  be  also  con- 
sidered as  unknown,  and  two  conditions,  kinematic  and  dynamic,  raust 
be  imposed  there.   For  simplicity,  however,  we  assvime  that  on  the 
free  surface,  (\/9)    is  small  compared  with  other  tenns  appearing 
in  (2),  and  that  this  surface  is  approximately  z  =  z   at  all  times. 
Since  the  pressure  there  is  p   ,  the  two  conditions  are  now  replaced 
by  a  single  linearized  condition,   (j),(x,y,  z  ,t)  =  0  ,  which  together 

TJ         O 

with  the  assumption  on   (\7<j>)    iraplies  that  (})(x,y,z  ,t)  =  constant, 
and  the  constant  may  be  taken  as  zero.   Thus 

(2.7)  <{)(x,y,z^,t)  =  0  ,    for  z  =  Zq  . 

This  approximation  was  discussed  by  Friedman  [1|],  and  is  commented 
upon  in  the  introduction  above. 

The  m.atheraatical  problem  vrhich  we  consider  is  that  of  finding 
ij)  and  F  satisfying  equations  (1)  through  (7),  given  F(x,y,z,0) 
and  P,(x,y,z,0)  (i.e.,  given  the  initial  surface  and  velocity  of 
the  bubble).   Once  ^     is  known,  the  pressure  field  can  be  computed 
from  equation  (2).   V/e  v;lll  assume  that  the  bubble  is  initially  a 
sphere  of  radius  A  and  that  its  initial  velocity  is  radial  and 
equal  to  A  ,  a  constant.   If  equations  (6)  and  (7)  are  omitted  (or 
If  z   and  H  become  infinite)  the  present  problem  reduces  to  that 
considered  in  [2],  provided  that  ^     is  required  to  be  regular  at 


i  . 
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infinity.   Hero,  as  in  the  case  of  the  reduced  problem,  the  assvunp- 
tion  of  very  special  initial  conditions  on  the  bubble  surface  is 
essential  for  applicability  of  our  method  of  solution. 

In  the  sequel,  it  will  be  convenient  to  employ  the  onorpiy 
balance  equation,  which  is  a  direct  consequence  of  the  preceding 
equations.   Since  in  the  linearized  problem  there  is  actual  mass 
flow  across  the  linearised  free  surface  (since,   ^  (x,y,z  ,t)  ^  0), 
it  is  not  at  all  obvious  that  the  total  energy  of  the  system  is  con- 
served.  This  is,  ncvcrthcloss,  true,  and  we  now  procood  to  shov;  it. 

Denote  by  B,  3^,  S^   the  surfaces  of  the  bubble,  of  the  free 
surface,  and  of  the  rigid  surface,  respectively,  by  V,  the  volwae 
of  the  water.   Let   dT  be  the  volume  element,   ds  -  the  surface 
element  measured  along  the  outward  normal.   Multiply  equation  (5) 
by  u  •  ds   and  integrate  over  B  ,  obtaining 


u  •  ds  =  0  . 


(2.8)  j  (Pli^t  ■"  |(V^)^]  -  Pq  +  P£Z  +  KV"^) 

B 
This  is,  as  we  shall  sho^^r  by  properly  identifying  various  terms,  the 

desired  energy  equation. 

Let   T  bo  the  kinetic  energy  of  water, 

(2.9)  T  =  |p  J  (V^^)^^i^  =  |p  J  ^V^  •  d G  =  i  p  j   ^y,|,    .  di  . 

V  S  B 

Hero   S  =  B  +  3^  +  S^  .   The  first  equality  follows  from  the  Green*  s 

theorem  for  potential  functions,  while  the  second  follows  from,  the 

fact  that  on   S-,  ,  (j)  =  0  ,  ;>:hile  on   S   ,  the  normal  components  of 

T/i   Is   i  =  ,|)  =  0  .   Now 
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where  V  is  the  velocity  ef  the  bubble  surface.  V/o  observe  that 

(V^  •  V^'t^  =  '^   •  (<l>tV4^)  -  k^A^   =  V  •  (^'tVM  .   Hence,  by 
Gauss'  Thoorera,  the  first  integral  above  is   2  |  '|)+.X7'{'  *  ^^   >   which 
rcdiices  to  2  j  <j>^\7'i>  *  ds  ,  since  on  3n,<{).  s  0,  v;hilc  on  S^  , 

"b 

h     =  0  .   As  to  the  second  term,  v/o  observe  that  on  the  bubble 
^n 

V  =  u  ,  so  that,  finally, 

(2.10)      p  ]l^^+   i(V*)^]^  *  "^"^  =  It  ^  =  i  P  It^  '  ^"^  *  '^'^^    ' 
B  B 

identifying  the  first  terra  in  (8). 

Let  ds  =  -  ds  ,  the  surface  cleracnt  measured  along  the  out- 
ward normal  to  the  bubble.   V/e  have,  since  u  =  v 

'         n    n 

!/ 
-^    -^     -^    -^'    d 
-  .  u  •  ds  =  j  V  •  ds   =  -tt:  V 

B 
Hence 

(2.11)  -  P^  )  ^ .  d  =  yvj) 

'B 
(2.12)  I^-Y  j  u  .  dt  =  -  KV-Y  1^  V  =  |^(4j  V^-"^)  . 


Lastly 

r  f  ^ 

(2.13)    +    zu  •  ds  =  -   !  zv  .  ds  =  -  ~   j  zdT  =  -  -J^  (BV)  , 

B  B  "v 

v/hore  B  is  the  location  of  the  center  of  gravity  of  the  bubble. 
Substituting  (10)  through  (13)  in  (8)  and  integrating  with  respect  to 
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timo,  wo  now  got 

(2,1k)  T  +  (P  -  pgB)V  +  -~  V-^"'^  =  constant  =  E^ 

O  V— X  o 


Tlao  energy  equation  (1[|.)  is  identical  in  form  with  that  ob- 
tained for  the  bubble  placed  in  infinite  water;  see  [2],  equation 
f20\  The  total  energy  E   is,  however,  not  necessarily  the  same  as 
in  case  of  infinite  water.   V/e  write 


(2.15)  E^  =  E  +  E   , 


where  E  is  the  energy  that  system  would  possess  if  the  water 
boundaries  wore  thrown  to  infinity,  while  the  bubble  has  the  same 
initial  data.   As  is  well  known. 


(2.16)  E  =  If   A^[|  pA^  +  P^  +  ^{^A^)-Y]   . 


III.  Moving  Coordinates  and  Dimensionless  Variables. 

The  method  of  solution  of  the  problem  posed  in  Section  II 
will  naturally  be  similar  to  that  used  previously  in  [2].   In  par- 
ticular, we  introduce  new  coordinates,  €,   >>,  ^,   as  previously,  by 

(3.1)  ?=x,   \=  7  ,        K=  z  -   B(t)  . 

B(t)   has  been  defined  in  Section  II  as  the  position  of  the  center 
of  gravity  of  the  bubble,  and  measures,  in  a  certain  sense,  the 


-8. 


migration  of  the  bubble.   The  moving  coordinate  system  defined  by 
(1)  is,  therefore,  such  that  its  origin  will  remain  at  the  center 
of  gravity  of  the  bubble  at  all  times.   V/e  also  introduce  the  polar 
coordinates  r,  9,  co,  relative  to  the  now  origin.   Since  the  problem 
still  presents  rotational  symmetry,  the  results  will  bo  indopcndcAt 
of  0)  ,  and  wo  can  assume  this  at  the  outset.   If  \<[g   write  the  equa- 
tion of  the  bubble  surface  in  the  form 

(3.2)     P(x,y,z,t)  =  r  -  R(0,t)  =  0  ,    or  r  =  R(e,t)   , 


equations  (2.I4.)  through  (2.7)  become,  using   *  =  ^ 

(3.3)     <j)^  -  -|  ^Q  =  R^  +  B(cos  9  +  -^i^  Rq)  ,    for  r  =  R(e,t) 

(3.i|)     Pq  -  p[g(rco3§+  B)  +  ^^  -   B(cos  9  ^^   -  ^^^  ^q) 


lf,2   ^  1   .2. 
2^K   ■"  I2  ^r' 


=  K(~   I  R^  sin  9  d9)"'^  ,    for  r  =  R(9,t) 


o 


(3.5)      cos  9  i   -  ^^^  ^   ^     =  0    ,        for  r  cos  9  =  -  (H  +  B) 


(3.6)     ({)  =  0   ,  for  r  cos  9  =  z  -  B  . 


In  addition  to  those     the  initial  conditions  on  the  bubble  are. 


(3.7)      R(9,0)  =  A  ,     R^(9,0)  =  A  ,    B(0)  =  0  ,    B(0)  =  0   , 


/4  '>. 


,,     -■  V  t'tiy'Sa.t.  r' 


■-_  I-.-  V 


aolcf 


.'lO.-u.^iJ; 


(f:,i) 


'"    V  <"•  ~) 
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while  the  condition  that  the  moving  origin  bo  at  the  center  of 
gravity  of  the  bubble  is 

(3.8)  \     R^   cos  e  sin  9  dG  =  0  . 

J 
o 


The  energy  equation  is,  explicitly. 


,2,„  ,     'r  ^  ,  sin  9, 


(3.9)       -  Tip   j  R  (R^  +  B[cos  9  +  -~~RQ])^(R,9,t)sln  9  d9 

o 

Id  It 

+  -^(Pq  -  pgB)  j  R^sln  9  d9  +  -^(-^  |  R^sln  9  d9)"'""*^=  E  +  E  . 


Finally,  wo  define  the  dlmensionloss  functions,  variables, 
and  constants,  r,  t,  9,  A,  b,  a  ,  a  ,  cT,  k,  k,  a,  hj   introducing  a 
unit  of  length  L  and  a  unit  of  time  T  .   These  quantities  are  de- 
fined as  previously,  by 

r  =  L  r  ,    t  =  T  t 

R(9,t)  =  L\(9,t),    A  =  La^,   A  =  LT'^a^,    B(t)=Lb(t) 
(3.10)/  |)(r,9,t)  =  LV^"^(?,9,t) 

&  =   LpgP"^  ,   k  =  '^-2. ,   k  =  max  k  =  J:-.(Xzl)r 

°  Y-1  ^         ^ 


[&  =  (^)V3    . 


For   the  linits  we   choose,  as  before. 


:."M 


Htn^ 


'•.     h'^il'l 


l!      "      '. 
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in addition  to  those,  we  define  two  new  dimensionless  ratios,   v 
and  [i   ,    by 


2 

(3.12)    V  =  Hpgp;^  =  E  ^  ^     ^  =   2^pgp;^  =  TT  "^    ' 


Clearly,  (J.  <  1  ,  and  equality  holds  only  when  the  pressure  at  the 
free  surface  vanishes.   The  ratios  v   and  [i     measure  the  distance 
of  the  bottom  and  top,  respectively,  from  z  =  0  where  the  bubble 
is  initially  located.   Thus,  for  fixed  cf  ,  the  larger  [x  and  v  , 
the  less  important  are  the  boundary  effects  compared  to  gravity 
effects. 

In  terms  of  all  these  new  variables,  the  previous  conditions 
become,  omitting  all  bars. 

Those  units  are  different  from  those  introduced  by  Taylor, 
Friedman  and  Shiffman,  and  frequently  used  in  literature  on  under- 
water bubble  motion.   If  we  denote  by  L'   and  T»   the  units  defined 
in  [I4.] ,    their  connection  with  the  present  is. 


-  =    /I     ^ 
2 


L  =  aL«  ,    T  =  /^     aT'  . 


The  reason  for  broalcing  with  the  custom  is  that  L  represents 
exactly  the  equilibrium  radius  of  the  bubble,  and  thus  has  a  greater 
physical  significance  that  L'   which  gives  only  an  upper  boumd  for 
the  maximum  radius  and  is  close  to  the  actual  maximum  radius  only  for 
high  bubble  energies.   The  new  unit  of  time  was  chosen  for  convenience, 
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(3.13)    A  k  =  0   y        for   r  >  R{0,t)   and  -(p  +  b)  <  r  cos  9  <  -^-b 

(3.11;)   g(9,t)  =   \>      -  \.  -  b  cos  G  -  -r^(:~  +  b  sin  0)  =  0  ,  r  =  \{9,t), 

(3.15)   h(9,t)  5  ^)^  +  |{[1  +  (x®)2](l|,^  +  i  sin  9)2 

+  2X^x^(^  ^g-^  ^   sin  9)  +  X^  -  b^j 

+  X   ii    \     ^^sin  9  d9)"^  +  cr(X  cos  9  +  b)  -  1  =  0  , 

k  '^  ^ 
o 

r  =  A(9,t)  , 


(3.16)  ^  =  cos  9  (})^  -  ^^g  ^  'I'q  =  0   ^   °n  r  cos  9  =  -  (^+  b)  , 

(3.17)  ^  =  0  .  ,   on  r  cos  9  =  -^  -  b  , 

f  h 

(3.18)  )  X^   cos  9  sin  9  d9  =  0   , 


whereas  the  dimensionless  energy  equation  is 

-ji: 
(3.19)    -  i     j  '^^i'^^   +  b[cos  9  +  ^^^  ^  XQ])^)(X,9,t)sin  9  d9 
o 

+  |(1  -  Ob)  j  \^sln  9  d9 
6 

+  a-^\(|  j  X^sin  9  d9)l-*^  =  ^-2(1  +  1)   . 
o 

The  functions  g  and  h  in  (ll].),  (15)  are  defined  by  those  equations. 
'Equation  (15)  is  obtained  from  (I4.)  by  using  (3)  to  eliminate  ^ 
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Our  problem  Is  now  to  find  (j),  \,  and  b  satisfying  (13-18), 
subject  to  the  initial  conditions 


(3.20)     X(9,0)  =  a   ,    XA9,0)    =  a^  ,    b(0)  =  b(0)  =  0  . 


a   and  a   cannot  be  prescribed  independently,  since  they  are 

related  by  equation  (2.16),  which  in  terras  of  the  new  variables 

becomes,  } 

I 

r 
t 

r 

(3.21)  [|  a^  +  1  +  k(aa^)"^Y]a3  =  ^'^   . 


IV.   Method  of  Solution. 

In  order  to  find  (j),  \,    b,  we  assume  that  they  can  be  expresses 
as  power  series  in  o"    (convergent  or  asymptotic)  in  the  following 
form: 


CO      00  /  j.n  \ 

(I^.l)   ^(r,9,t)  =T:<^     11  oit)ir-^'''-'^h(cos   9)  +  4);:(r,9)) 


CO      00 

£  ^   E 

n=o   m=o 

CO        CO 

Z   ^'Z 

n=o    m=o 


-Z   ^T  c_(t)j 


nra 


CO     ^  CO 


(i;.2)   x(9,t)=j;;;  cf^x;  a   (t)p  (cos  9) 


„_^    ^ imi'  '  m 

n— o    m— o 


(1^.3)   b(t)  =  2_.  b^(t) 


CO 

n=o 


Here,   Cj^(*)  >  ^nm^^^^  ^n^^^   ^^®  flections  of  time  to  bo  determined. 
In  view  of  (3.20),  they  must  satisfy  the  following  initial  conditions: 


•13- 


C 

'a  (0)  =  a  ,   a  (0)  =  a 
oo       o     oo       o 


(U.W    -j  a^(0)  =  a^^(O)  = 


=  0   for  n  +  ra  >  1 


1  b(0)   =b(0)=0    for  all  n   . 
*y^  n       n 


We  require  that  the  \"     be  regular  harmonic  ftmctions  irx  the 
5trip  -(^+b)<^<-^-b,  such  that  J^  =  r"^"^'^-'-^P^(cos  0) 
+  \'^     satisfies: 

f  1^=  0  ,    for  ^  =^-  b 


^v\ 


(i;.5) 


& 


?^  =  0  ,    for  ^  =  -  (^+  b)  . 


Then  (1),  which  certainly  satisfies  (3 •13)  will  also  satisfy  tho 
boundary  conditions  (3.1<3,  3«17),  no  matter  what  the  c    turn  out 


to  bo. 


I  ""  I  ■■"■■ 

The  ^       are  determined  in  the  Appendix  A  .   Obviously  <p" 


depends  on  (f    and  t   (through  b)  sincoj  as  equations  (5)  shox^;, 

their  boundary  values  depend  on  o'  and  t  .   Since  ifc'"'  is  regular 

"^m      ^' 

at  r  =  0  and  has  rotational  symmetry,  it  can  be  represented  as  a 
linear  superposition  of  axially  syrametric  spherical  harmonics  regular 
at  r  =  0  ;   hence 

k— o 

The  coefficients  d,  ,  (see  Appendix  B),  depend  on  cf    and  t  ,  and 
as  is  shown  in  Appendix  C  ,  they  can  be  expanded  in  power  series  in  (f 


-11^. 


00  // 

On  the  other  hand,  since  the  problem  reduces  to  that  of  the 
spherical  bubble  theory  when  o^=  0,  we  have 


ik'Q)  a^_  =  c   =  b  =  0  ,    for  m  >  0  , 

^  om    om    o     *  ' 


whereas 


ik'S) 


a^^(t)  =  a(t)  >  0 
c^^(t)  =  -  a^a  , 


where  a(t)   is  the  well  known  zero  order  solution,  discussed,  in 
particular,  in  [2]. 

In  view  of  (7),  (B)  and  (9),  ^,    <j)^,  \>q,    ^^     for  r  =  X(9,t) 
can  be  expressed  as  power  series  in  a" .   Likewise  the  volume  of 
the  bubble  can  be  expressed  as  power  series  In  cT"  .   To  find  the 
b^,  a^^  and  c^^^  for  n  >  1  we  substitute  the  power  series  thus 
obtained  in  (3.1'-!-)  and  (3.15),  and  equate  the  coefficients  of 
crp^(cos  9),  for  every  n  and  m  .   This  way  we  get  a  sequence  of 
equations  Involving  the  coefficients  and  their  derivatives. 

To  obtain  these  equations  explicitly,  it  is  convenient,  as  was 
done  previously  in  [2],  to  introduce  the  functions  g(9,t)   and 
h(9,t)   defined  by 

_  00        .         00    .        * 

(1^.10)   6(9,t)  =  6(9,t)  -f  g^  a-^(b^P,  +  E^Ca^^f  2|  a^ 

+  (m+i)a"^^^2^c  ]p  )  5  r;  cj^  r:  g  p 

nm-"  m'       j—  ^     ^nm  m 

n— o    m— o 


-.'■■•qa 


S  O     -i  ■ 


ti:  ■ 


-15- 


00  ,  on                      CO       /     d     \ 

(il.ll)      h(e,t)   =  h(9,t)    -Z:     <:^  (-3y^  ^^^^\o  ^II   "^^2        ^nm 

n=l  k  m=o      a 

/      ,  -,  X  .  CO           ^     CO 

+aa       +   a-^-'^'-^^c  ]P    )  =  V"     c/^  V     h     P      . 


Now,    inserting    (1)    through   (3)    In   O.lq.),    we   obtain,    for     n  >  1 


,    c        =   -  -^T  a"'"*'^ta       +   2-2  a        -   g      ]    ,        m  =1=   1 
ran  rn+1  nni  a     nm       '^xm.     '  ' 


CT=--ia-^[b     +aT+2-aT-g,],     m=l. 
nl  2  n  nl  a     nl        ^nl 


Similarly,  inserting  (1)  through  (3)  In  (3.15)  and  making  use  of  (12) 
we  get,  after  some  manipulations, 

^^no  ■"  2^^no  +   (^  -^  3yx  ^~'^~    )^^o  =  ^^^o  ^  ^ag^^  +  h^^   ,   m  = 
(1+.13)     ^    (a\)*    =    (a^tg^l   -   a^iJ)*   +   2a\^   ,        m  =  1 

•  •  ••  ••  •  • 

aa       +  3aa       +    (l-rrjaa       =  ag       +    (ra+2)ag       +    (m+l)h     ,    m  >  2 
nm  nm        ^  nra         *^nm        ^         '    '^nm        ^         '   nm'        — 

So  far,   g  and  h  are  general  formal  power  series  in  a'    with 
coefficients  depending  on   a^^,  c^^,  and  b^  .   If  g^^  and  h^^ 

were  independent  of  a    for  m  +  1  ,  and  of  b   for  m.  =  1  , 

■^  nm         '  n 

equations  (13)  together  with  the  initial  conditions  (Ij.)  would  deter- 
mine the  a    and  b   \iniquely,  provided  a  ,   wore  also  knoi^fn. 
nm       n        ^ ,   ±  j^j^ 

Even  so,  in  general  we  would  have  to  solve  for  each  n  an  infinite 
set  of  equations,  and  therefore  we  would  not  be  able  to  get  beyond 
the  terms  of  first  order  in  cr'  .   Fortunately,  the  following  theorem 
can  be  proved: 


i:u.)    i-r^ 


-16- 


1.  for  each  n  ,   g   =  h   s  0   for  m  >  n  . 

2.  for  each  n  ,  the   g,^,  h^^  depend  only  on  the   a^^,  b^,,  Cj^ 
with  k  <  n 

A  similar  theorem  was  proved  in  [2],  and  since  the  method  of 
proof  is  the  same,  wo  omit  it  here.   (The  theorem  proved  in  [2]  was 
somewhat  stronger  in  that  we  managed  to  show  the  g   s  h   =0 
also  whenever  n  +  m  =  odd  integer.   Tliis  is  not  true  In   the  present 
case. ) 

Prom  the  first  part  of  the  theorem  it  follov:s  then  that  the 
equations  (13)  are  linear  and  homogeneous  for  m  >  n  ,  and,  since 
all  the  \inknowns  are  subject  to  homogeneous  initial  conditions,  that 


(I4..I/4.)  a^,^  s  0   for  m  >  n  . 


Hence  for  each  n  ,  wo  have  onl^y  a  finite  nurabor  (exactly 
n+1)  of  equations  to  solve.   Also  using  (ll|.)  in  (12),  one  gets 


(I1.15)  c   SO   for  m  >  n  . 

^    '  nm 


Thus,  the  expansions  (1)  and  (3)  simplify  to 


CO     n_  r  4-1  ■» 

('4. 16)  ^(r,P,t)=2__  c/"2_     c   (t)(r-^^'"-'^P  (cos  0)  +  i)^) 

n=o    m=o 


CO     n 
n=o    m=o 


ik-n)  7i(G,t)=^  </*  T;  ^nm^^^^m^''°^  ®^  • 


In  view  of  tho  second  part  of  our  theorem,  the  equations  (13)' 


-17- 


for  a  fixod  n  and  n  <  n  ,   n  ^  1  ,  can  be  solved  provided  that 
we  know  a,  ,  b,,  c^^  ,  with  k  <  n  .  Also,  since  by  (3.19)   a  ^^ 
can  be  detorr;;incd  In  torris  of   a^^  with  k  <  m  ,  the  equation  for 
b   can  be  likewise  solved.   Hence,  the  systora  (13)  can  bo  solved 
in  succession. 

V.  Bubble  Coefficients  up  to  Temis  in  d-^    . 

Up  to  terms  in  d       the  expansions  for  (|),  X,  and  b  arc 

(5.1)  ^  =   c^^(r-l  f  d^^  +  d^3_rP^  +  d^2^%)  +  [c^^(r-l+  d^^+  d^^rP^) 

+  c^i(r"%  +  ^lo)]^+  ^^2o(^"^-^  ^00^  ■"   ^21^"%  ^  ""zz'^'^^Z^' 
•*■  ^^30^"^  ■"   "31^'^^  "-   "32^"^P2  -^  '^33^'S^^^ 

(5.2)  X  =  a  +  o.^^cr  +   (a^^  +  !X^^?^)d^  +   (a^^  +  a22P2  "^  a  P^)cr^  +  ... 

{S'3)      b  =  bj^  o'  +  h  d^  +   b  d'^   . 

In  the  expansion  for  X     tenns  Involving  a^^^  wcro  omitted  since, 
as  follows  from  (3.18)  they  are  zero  for  n  <  k   >      (This  assertion 
is  not  true  for  all  n  ,  though.   As  may  be  verified,   a^, 
-  -  Yq  a  °-22°'33  T  0.) 

In  this  section  wo  consider  the  equations  for  a    with 

ITTl 

m  ^  0.   The  equations  for  bubble  coefficients  with  m  =  0   are 

derived  in  Section  VI  in  a  simpler  way.   Since,  however,  the  present 

equations  will  depend  on  the  a^^  ,  some  of  the  results  of  Section 

VI,  naracly  those  for  a^   and  a^  ,  are  already  used  here. 


-18- 


Lot  a  ,  b   denote  the  bubble  coefficients  when  the  boiin- 

nr;r   n 

darles  arc  rcrr.ovod.   The  equations  for  these  were  found  in  [2],  and 
will  be  omitted  here.   Let 


a   =  a   +  a 

nil    nm    nm 


(5.i^) 


b   =  b  +  b 
n     n    n 


Thus  the  a    and  b   are  corrections  to  the  bubble  coefficients 
nm       n 

due  to  existence  of  boundaries.   It  was  previously  established  that 

For  n  =  1  ,  the  only  coefficient  (except  for  a,  )  to  be 
determined  is  b   .   V7e  find,   g-, -1  =  0  ,   h, ,  =  a  ,  and  consequently 


(5.5) 


b^  =  0  , 


'11 


1  3* 


Hence  b^  =  b,  .   It  was  previously  shown  that  the  formula  for  b. 

Is  equivalent  to  the  Herring  rise  formula. 

For  n  =  2  ,  the  only  coefficients  (except  for  a^  )  to  be 

-  ^  x-  2o 

determined  are     bp   ,    ajtid     a.p^   .      V/e  find, 

^  621  =   -    ^^\   ^"""^10  ^   ^^^  ^001^ 


(5.6) 


S22  =  0 

^^21  "  ^^10   -   ^(^^^)    ^ool 


V  ^22 


3  12 

11  \        ' 


and,  on  substituting  in  ([j.l3),  obtain 


i^v/:) 


-19. 


(5.7)         ^2  =  ^2  =  ^2^001  ^  ^2^012  ' 


where   P„   and  p„   arc  dctorminod  by  , 

j  (a^P^)'  =  I  b^a(ia^+  [aa  -  i^]  f  a('u)d'r)  ,  ^^{0)=h2(0)=0 


(5.8) 


(a^P2)'  =  -  3a^(a^i)*  ,   P^CO)  =  p2(0)  =  0   , 


while. 


(5.9) 


^22  -  ^  * 


In  terms  of  bp  and  ^ppf   ^"^o  S^t  from  (6)  and  (i|..12)  , 
1  °21  =  -  |f^"^^2  ■"  ^^l^^^lo  ■"  ^^^  ^012^ 

(5.10) 


I      _    1   2,  2    X* 
Cpp  —  -  TT  a  ^a  ^pp  J    • 

In  establishing  equation  (8),  v;o  used  the  result  a,   =  a-, 
=  -  ~  d  -a  [  a(T)dT  ,  (soo  equation  (6.12)  )  . 

Per  n  =  3>   }    the  only  coefficients  (except  for  a^  )  to  be 
determined  are  b-,  a-Do*  and  a,,^  .   We  find 


r 


(5.11) 


'31 


32 


a-^[3ft2^2^^^  +  3l>iaa2^  +  3b^a2(a2^  -  ^a^^)Ha\^^)d 


ol 


-   a 


•^ro/.3 


[2(a^a^^a22)*    +  ^^^^lo'^'zz  "^   2a'^^^o23^ 


^33 


h 


31 


18  ;    -1 

"T  ^^      ^22 
3     -2//      2 


1   2 


=  I  a  ^(b-j_a^a22)   +  3a2Q  +  3a  '  aj^-    [(a-^)    a^^+   a(a^a^^)]d 


ol.- 


-20- 


^32  =  -  2  ^^V  ^^^^ol2^  -  ^^(^^^^*^o23  ■"  ^lo'W^^lo^  i;a-^aa^^)a22 


+  (a^^^  +  i+a"  aaT^  +  2a"  aa^^ja 


^lo 


lo'"22 


3   *   -1    '   9 
^33  =  -  5  ab^(a  a^^)     ^  -^  ^^Z   ' 

Substitution  in  (i|..13)  now  yields,  using  a,  ,  a   ,  dGfinod  by 

a,   =  a,  d  ,  ,   a^  =  a.„  d  T,(soe  equations  (6,12)  and  (6.11^.), 
lo    lo  ool  *    2o    2o  ool' ^     ^  \        T,j 


(5.12) 
(5.13) 
(5.11^) 


b-  =  p-d^  ,  +  P-d  ,d  ,^ 
3    3  ool   "^3  ool  ol2 

/v  t         »» 

32    32    iZ   ool    32  0I2    32  023 

a^^  =  0  , 


t         t         " 

where  p   ,  p   ,  <^^oy    "^to  *  ^^^  ^-op     satisfy  the  following  equations: 


(5.15) 


(a^p^)'   =   -   3(p2a^a^^  +  ^laa^^  +  b^a^a^^)* 

+   6(a^a2Q  +   aaj^)    ,      p^(0)   =  p^(0)   =  0 

(a^p')*   =  -  3(p2a^a^o''   ^^^o^'^' 


•  I 


+  12a^a(a^a^^)    ,        P^(0)   =  p^(0)   =  0 


•  • 


(5.16) 


r  •  •  ••  ••  O**  ** 

aa^2  "^  3aa^2-   ^^^32  "  "   2  ^1*^2  "    ^^^10^22  "^  ^a^  ^a^^~   a,  ^a. 


lo"22        lo"22 


\        •  .t  •  •  t  •  •  t 

1     aa,2+  3aa^2  "   aa  2 


a^2(0)   =  a^2(°)   =  ^ 


2  ^1^^2  "*■  ^^^^    '      a^2(°^   =  a^2^°^" 


aa^2+  3aa  2  -    aa       =  -   5(a*'a)      ,      a      (0)   =  a.-,^{<:))   -  0    . 


V 
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Forraulae  for     c-,,    c-„,    aiid     c^-,      can  bo   obtained  in  terns   of     b., 

31   32'       3j>  3-L 

a^P  and  a-,^  using  (6)  and  (4. 12).   Howovcr,   c^-  and  c^^  are 
of  no  interest,  and  the  expression  for  c-^  is 

(5.17)   C33_  =  -  |[a3b3  +  3i^a\^   +  3b^aa2^  +  ^K^^^^Zo   "  j   ^22' 

^  ^^\o^*^ol2^   • 


For  the  equilibrium  bubble   (a  =  1)   discussed  at  length  in 
[2],  wo  find  a,   =  0  ,  a^  r  0  ,  and  consequently  all  a„_  ,  b^ 
equal  identically  to  zero  except  for  a-,   .   Thus,  tho  only  case 
which  can  be  solved  explicitly  is  uninteresting  from  the  point  of 
viev;  of  effect  of  boundaries.   To  obtain  quantitative  results  in 
the  general  case,  equations  (3),  (15)  and  (16)  must  be  integrated 
numerically. 

VI.   Use  of  the  Energy  Equation  -  Bubble  Coefficients  with  m  =  0 

The  bubble  coefficients   a    may  be  computed  by  solving 
equations  ([|..13)  with  m  =  0  ,  namely 

(6.1)    a'4^  +  3a;^,  +  (*a  f  3t|  a'^^Y+Dj^^^  =  ^'^^^  ^  sig^^  ^  ^^^   , 
Using  the  equation  defining  a(t)  ,  namely 


(6.2)         I  a3;2  +  a^  *  ,r^—   a'^lv-D  =  £-3 

^  k(Y-l) 


fi'  I' 


.     ) 


-22- 


(cf.  [2]  or  earlier  works  on  spherical  bubbles),  one  verifies  that 

• 

a{t)  is  a  particular  solution  of  the  homogeneous  equation  associated 
with  (1).   Formally  then,  the  problem  of  solving  (1)  may  be  reduced 
to  two  quadratures,  since  there  exists  a  formula  involving  two 
integrations  which  will  produce  the  general  solution  of  a  second 
order  linear  differential  equation  once  a  particular  solution  of 
the  associated  homogeneous  equation  is  known.  As  is  shcxm  In 
[5]>  this  formula  involves,  in  the  case  when  this  particular  solu- 
tion vanishes  at  some  points  -  which,  indeed,  is  the  case  with 
a(t)  -  the  computation  of  finite  parts  of  divergent  integrals, 
unless,  of  coiirse,  the  integrations  can  be  carried  out  explicitly. 
V/e  shall  show  here  that  _in  pai^t  these  integrations  can  be  carried 
out,  but  in  general  it  is  preferable  to  solve  (1)  numerically  rather 
than  to  use  formulae. 

In  any  case,  equation  (1)  can  be  reduced  to  a  first  order 
equation  which  is. 


(6.3)    kk^^.   aa^^^  =  k^^  =  a'^  j  (ag^^^  +  zkg^^   -f  h^^)a2;dt,  a^_^(0)=C 


To  proceed  further  we  must,  of  course,  evaluate  g    and  h    in 
^  *=no        no 

terms  of  the   a,  //  ,  b   v;ith  k  <  n  ,  and  the  actual  expressions 

for  these  are  a  great  deal  more  involved  than  the  formulae  for 

g    and  h    with  m  >  0   encountered  so  far  (equations  5.6,  11). 
^nm        nra  \  m        >   >    / 

A  simplification  is,  fortunately,  possible.   In  Section  II 
we  derived  the  energy  equation  which  in  dimensionless  form  is 
equation  (3.19).   The  energy  equation  is  a  first  integral  of  the 


problem.   Hence  on  substitution  in  it  of  our  expansions  for  ^,  X, 
and  b  ,  and  on  identification  of  coefficients  of  various  powers  of 
(f  ,    one  ought  to  get  a  sequence  of  first  order  differential  equa- 
tions for  the  bubble  coefficients.  Now,  since  the  equations  for 

a    with  m  >  0  cannot  be  reduced  to  first  order  equations,  we 

nm  ^ 

expect  to  get  from  the  energy  equation  the  equations  for  a   , 
namely  the  equations  (3).  Moreover,  since  the  integrations  intended 
in  (3.19)  are  with  respect  to  Q  and  not  with  respect  to   t  ,  and 
since  the  expansions  used  are  algebraic  forms  in  the  bubble  coeffi- 
cients and  their  derivatives,  the  equations  so  obtained  will  involve 
the  bubble  coefficients  algebraically.   Thus  we  shall  not  only  avoid 

the  computation  of  the  g   ,  h   ,  biit  will  also  obtain  the  integrals 

no   no  ° 

k    in  (3)        ■  ... .    g^3  algebraic  exor  ess  ions  in  the  bubble  co- 

efficients.   (Prom  the  form.  (3)  of  k    it  is  of  course  not  at  all 

no 

obvious  that  these  integrals  can  be  evaluated  explicitly. ) 

We  now  proceed  to  express  equation  (3.19)  in  terms  of  bubble 
coefficients  up  to  order  I4.  in  o'  .   Let 


(6.I4.)         a„  =  E  <<"a.^  . 

n=o 

In  view  of  the  theorem,  of  Section  IV,   a  =  O(o^)  ,   c  =  0((/'^)  .   In 
place  of  a   it  is  convenient  to  uso  the  average  radius  of  the 


00 

c„  =  T" 

m       ^— 

nm 

n=o 

bubble,   a  , 


(6.5)  a  =  (|  J  X^  sin  9  dQ)-^^^ 


in  terras  of  which  a   is,  up  to  terms  in  d^   , 


^/i;,^- 


f  i 


ir         t!u 


\  ■. 
(  -■ 
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,LL\  1-12 

(6.6)  a=a--paa_. 

o        p      ^ 


Using  the  above  notation,  and  substituting  (1^.16,17)  in  (3 •19),  we 
now  get,  on  carrying  out  tho  intended  integrations, 

■3       'l*         P*  2'1      -1    * 

(6.7)   -^rc(aa--^  a^a^  +  d^^a^a)  +  c^d^^a  a.  +  -^  o^{a     c^) 

+  b(ic  d  ,a^  +  ic,  -  |c,a"-a.J]  +  (1  -  crb)a^  +  ka'^Ya-3(Y--"'-) 
30  ol     3-'-   j?l    '^ 

=  a"3(l  +  I)   . 


V/e  still  need  to  express  c  ,  c  ,   and  c-  in  terras  of  bubble 
coefficients,  of  x/nich  c   is  needed  to  order  l\.,    c-,   to  order  3  > 
and  c    to  order   2   in  cr  ,    if  (7)  is  to  be  correct  to  order  k 
in  c'  .   c   can  be  evalueated  exactly  (to  any  order)  by  evaluating 
the  surface  integral  of  the  normal  derivative  of  ^     over  the  sur- 
face of  the  bubble.   Since   (j)"   (see  equation  (A.9))  is  a  regular 
potential  function  inside  the  bubble,  its  contribution  to  the 
surface  integral  is  zero,  b;/  a  well  known  theorem.   Hence 


Bubble       Bubble   "'^ 

7t 


=  -  2:c  >_  c^^  I  [(ra+DP^^^  +  ^''^\-^l^y'^\   sin  9  d9  . 


m 


The  integrand  above  is  -  'Ig^''^"  ^n  ^^'^   ®^  '  ^^     ™  +  0  >  v;hile  for 
m  =  0  it  is  eqixal  to  sri  0.  On  integrating  one  gets  then. 


q  =  -  i|7tc^ 
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On  tho  other  hand,   (^  can  be  computod  by  using  the  bovmdary  condi- 
tion at  the  bubble,  equation  (3.1^!),  natnol:/ 

Q  =  27t        [A^   +  b(coa   0   +  -^^.   sin  9)]X^sia  9   d©  ^- 

71 

TT  / 

=  2it     ('    ~(X^)*sin  e  +  •?:{)  4o(>^'^sin'^0)d9  =  -^  -|r(-i       X^sin  QdG) 
I     j>  2     c'.y  _5   a.!;  d     i 

o  o 

2* 


^equating  the  two  expressions  for  Q;  ,  one  gets 


2* 
(6.3 )  c  =  -  a  a  . 

^    '  o 


The   expressions   for      c,      aiid     c^      to   the   desired  order   can  be   ob- 
tained by   conbining  rooults    (5«5;10).      One  gets    thus. 


T'OOO  1'-*                                                                                                    1 

(6.9)  c     =   -  -  b(a-^   -  -^  a^a^)  -  -  a-^ad  -,    ,      to   tonns    in     cr 

1      2      2         *  ? 

(6.10)  c„  =   -  -:r  a    (a  a^)  ,      to   terms   in     c^ " 

t-  5  <- 


Substituting  (8),  (9)  and  (10),  in  (7)  ,  wc  get  finally 


(6.11)  I  {;2(a3  +  ,.J-^   +  id^^d.^a^  -  f^aa^j^  b2(la3  -  1^\^) 
-f  (1  -  c-o)a3  -.  i,^--3r^-3(Y-l)  =  --3(1  +  I 
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Tho  last  part  of  this  equality  was  obtained  by  evaluating  the  energy 
at   t  =  0  and  using  (3.21).   The  equations  for  the   a  ^  are  now 
obtained  by  replacing  a,  a„,  b  by  their  expansions  in  (i^)  and 
identifying  the  coefficients  of  various  powers  of  &  .      Results  arc 
as  follows. 

1.  For  a^Q(t)  =  a(t)  we  got  (as  expected)  equation  (2) 

2.  For   a,    wo  get  the  equation 

^^10  -  ^^10  -  2  ^ooi(-~r  -  ^^ )  • 

a 

•  • 

This  equation  can  be  solved  explicitly  if  a  =  0  ,  while  if  a  =f=  0  , 

the  solution  is  obtained  in  terns  of  finite  parts  of  divergent 
integrals.   V/e  shall  asstxtne,  hero  and  subsequently,  that  a  =  0  .'" 
One  gets  then. 


(6.12)  a^^  =   a3_^d^^3_  ,     c^^^  =  -  |  a  j  adT 


3.    For  ap   we  get  the  equation 


1              1*2        la2  1*2          *' 

aa^^   -   aa„^  -  rr  h^    -  rp^  b     +  —  -r  a.,  -  -^  a,      -   afva,  ^d     , 

2o             2o        3      1        12     1        2          lo  2     lo               lo   ool 

n. 

1  *2 

2  lo   ool 


The   significance   of   this   assrunption  was   discussed  in    [2}    • 
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Using  the  formula  (12)  for   a,   ,  one  now  finds  that 


2o    2o    2o  ool  ' 


where   a^^   is  the  part  of   a^   found  in  connection  with  the  bubble 
do                     ^  2o 

in  infinite  water,  (sec  equation  (5J^)  in  [2]),  Xirhile  a^   is  given 
by 

t  t  t 

{6. Ik)        a^Q  =  -^  {a(  I  adT)  +  2aa  j  adi  +   a  j  a  dt  }   . 


k<  For  a-   we  get  the  equation 


^^30  -  ^^30  =  a  ^lo^2o  ^  ^J-R   ^    ^lo  "  ^lo^2o 

-  I  ^(^10  -^  ^ii^o^^ool 

•  .  2.'  P  1  *  * 

-  aa,  ^a,  ^d^^-  -  -a  (2ab,  d^,„  +  a^  d  ,)  -  -7  b-b_ 
lo  lo  ool   2      1  0I2    2o  ool    6   12 

1  2*  *        1 

-  2  ^  ^^1^012  ■'3^2* 


Using  previous  results  for  a-,  ,  ap  ,  and  b^  ,  one  now  finds  that 

(6.15)        a-,  =  a_  d  t  +  a_  d  t  „  +  a_  d-^  ,  , 
^    ^'         30    30  ool    30  0I2    30  ool  * 


t 

where  a_   and  a_   satisfy  the  equations 
30        30        <}  ^ 

t  t 

(6.16)   aa^^  -  aa^^  =  ■^(a  j  ada  -  aa)a2^-  -^(a  /  ada  +  a  )a2 

o  b 
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(6.17)  aa^^  -  la^^  =  -(a^ab^  +  aa^b^  +  ^  h^li'^   "3^2^   * 

with  the  Initial  conditions  a^  (0)  =  0  ,  a,  (0)  =  0  .   (For  the 

30  30 

definition  of   p   and  (3   see  equations  (5.3).)   The  equation 

ri 

that  a,.   satisfies  can  be  integrated  explicitly  and  one  gets 

t  t      t 

(6.18)  a^^  =  -  -i^  {  ^*a(  ;  ad-r)^  +  I(  |  ad'c)(  (  a^d-c) 

o  6      6 

t  t        t        t 

+  (2aa  +  a  )(  j  adT)   +  3a  a  '  adT  +  aa  ;  a  dT  +  a  i  a-^dT}  . 

-''  i  J  J 

o  000 

The  equations  (16,  17)  arc  singular,  since  the  coefficient   a  of 

the  leading  teira  vanishes  at  some  points.   Yet,  as  we  Icnow  already, 

they  can  be  solved,  although  it  is  not  convenient  to  use  then  in  the 

given  form  for  n\americal  solution.   We  know,  however,  that  the  a 
•^  '        '  no 

satisfy  the  equations  (It.  13)  v;ith  m  =  0  ,  which  are  regular,  and 

therefore  a^   and  a.   satisfy  equations  T>rith  the  same  principal 

part  and  finite  non-honogeneous  terra.   Actually,  equations  (lj..l3) 

can  be  obtained  from  (16,  17)  by  multiplying  them  by  aP   ,  differen- 

2' 
tiating,  and  then  dividing  by  a  a  .   One  gets,  however,  simpler 

regular  second  order  equations  by  differentiation  of  (16,  17)  vjithout 

3  • 

prior  multiplication  by  a   .   Tlie  result  is,  on  dividing  by  a  , 

(6.16)'   a3^  -  (f)a3^  =  -  |{(^)*  (  ad-u  ^  2a(^)  +  ki}\^  -  ]  \ 

a  a  J         a 

t 

+  I   b^(a'^a  (adT  -  |)  +  a'^b^p^'  ^=^30^°^  "  ^30^°^  "  ° 
o 


I 
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t**t  flt  Ofp  ••  p  **l 

(6.17)      a        -    (t)^.^  =  -    (2aa  +   a   )b^  +   aab      -    a""  -   ab^p      , 
a     -^ 

The  right  hand  terms   in   (16)      and   (17)      were   arranged  so   as   to   ex- 
hibit  their  finiteness.      The   only  doubtful  parts   are   those  involving 
•••  • •• 

a       a 

—  and  (~  )  .   Prcm.  equation  (2)  one  finds, 

a       a 

^  =  I,a-\-^(l  -  f(3Y+l)k[aa]-3^^-l)) 
a  ^' 

•  •• 

(^  )*  =  -  20r\-^(l  -  J^(3r+l)(3Y+2)k[aa]-3(Y-l));  , 
a 


which  are  both  finite  for  all  t  sinco  a  >  0  . 

VII.     Discussion  of  Results. 

Formulae  for  the  bubble  and  migration  coefficients  are  too 
involved  to  allovi  for  a  coiapletc  discussion  of  our  results.   In  the 
companion  report.  Underwater  Explosion  Bubbles  IV,  we  use  these 
formulae  for  computations  in  a  special  case.   Here  we  shall  carry- 
out  only  a  qualitative  discussion.   In  particular,  we  shall  show 
that  certain  effects  predicted  by  other  authors  are  easily  deducible 
from  the  present  study. 

V/e  recall  that  to  the  third  order  in  o^ 

fx  =  a  +  a^^o^+  (ag^  +  ^22^2^^  "^  (^3o  "^  ^32^2  "^  ^33^^'^'^ 


I  (7.1) 


where 


") 


b  =  b^0'+  b  cf^  +  b^cT^ 
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/■ 


(7.2) 


lo    lo  ool 


2o    2o    2o  ool 


^22   ^22 

^30  =  ^30^001  ■*■  ^^30^012  ^  ^3o%ol 

^32  =  ^32^001  ■"  ^32^012  "■  ^32^^023 


a^^   a33 

^2  =  ^2^001  ^  ^2^012 


b_  =  b_  +  p^d'^  ,  +  P-d  ^d   T^   . 
3    3    3  ool    3  ool  ol2 


The  barred  quantities  are  those  that  would  be  obtained  if  the 
boundaries  were  removed  to  infinity. 
1.    Theory  of  bubble  periods. 

In  vlevj  of  formula  (6.6),  the  average  radius  of  the  bubble  a 
is  to  the  third  order  in  cT"  given  by  a  , 


(7.3) 


a 


a_  =  a  +  a.-,  .o^+  a^_(5^  +  a.~   6" 


Uo 


20 


^3o 


In  general  a  is  not  a  periodic  fiinction,  but  may  be  assumed 
approxiraately  periodic  for  small   o^ ,  since   a(t)   is  periodic.   Let 
'U  be  the  half  period  of  a(t)  .   The  half  period  t(o')   of  a  can 
be  defined  as  the  time  elapsed  until  the  first  maximuin  is  reached. 
That  is,  i{cf)      is  defined  by  the  equation 


-31- 
(7"^)  i(T(cr))  =  0  . 

Assume  that 


(7.5) 


T(d')  =  T  +  T^o^+  ...   . 


substituting  in  (k) ,    expanding  and  keeping  terns  up  to  the  first 
order,  we  get 

^^•^^  ^("^^-^  o-CT^'aV^-)  +  ;^^(^))  + ...  =0   ^ 

Solving  for  -u^^  ,  and  using  (6.12),  we  get 

''-'^  '^l=Kol(/a(-)^-)=K,,(  r^^-fda)   . 

o  Q^  a. 

"Tnin 

This  formula  is  Identical  with  that  obtained  in  [1,],  p.  i8i,  fernula 
(1.15). 

The  formula  for  d^^^  ^hich  is  derived  in. Appendix  C  is 

^^•^^         ^001  =  -  J^i^M   +  log  2)  , 
where 


(7.9)         f(x)=2xf-  LiDl  .  -  v-LL 

s=o  (2s+l)'^-x'^  v+M- 

The  function  f(x)  has  been  tabulated  in  [^4],  p.  190.      ^ 

ool   """^ 


( ■•', 


'  v'..  <. 
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negative  for  x  >-  ^  ,  and  positive  for  x  <  -  -r   .      Thus  the  bubble 
period  is  smaller  or  larger  than  that  of  a  corresponding  bubble  in 
infinite  water,  depending  on  whether  n  <  2v   or  |J,  >  2v  . 

2 .     Theory  of  jnip;ration . 

To  the  first  order,  the  migration  of  the  bubble  is  given  by 
the  Herring  fomula, 

t         T 

(7.10)  b  =  crt),  =  20^  (  -—   f  a^(o')dcf  . 

^      J     a^(T)  J 

For  derivation,  see,  e.g.  [2].   Experimental  data  indicate  that 
Herring's  formula  predicts  a  too  large  migration,  and  the  discrepancy 
is  of  the  order  of  50°/o  for  moderate  explosions  (300  lb.  TNT), 
For  bubbles  in  infinite  water,   bp  5  0  ,  and  b-  <  0  •  This  t^hows  that 
our  results  are  in  qualitative  agreement  vrith  experiraental  data. 

In  presence  of  boundaries,   b„  -^  0,  and  it  dominates  the  third 
order  corrootion.   We  shall  estimate  its  effect  hj   computing  b_  at 

« 

the  time  of  the  first  minimum.   In  this  estimate,  we  assume  that  b, 
is  appreciable  only  during  a  short  time  interval  near  the  first  mini- 
mum of  the  bubble,  and  that  the  bubble  radius, is  during  most  of  its 
m.otion.near  its  maximura.   Then  formulae  (5«8)  yield,  approximately, 
at  t  =  2t 

m 

(7.11)  ]    ., 
p^  ^  12a;^^'T:A 

Here,  M  refers  to  data  at  the  maximura,  m  ,   to  data  at  the  mini- 

k'2 
m.um  of  the  bubble,  and  A  is  the  average  value  of  a^a  .   Both 


IX 
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•      •  t 

Pp  and  p^  are  positive.   Using  the  above  estimate,  we  noi-j   get 

(7.12)1  MSt)  '-  h,cr  +  ^^d^  '-^  l^{Jl)^o'i  +   3(~)b,  (2T)a  6^id     , 
J  X  d  a^        a^  1     ra    ooi 

+  12  a'^  Acf^Td^T, 

m       012 


For  X  >  -  TT  ,    i.e.  for  [x  <  2v  ,  both  d  ,   and  d  ^^     are 
J  ool        ol2 

negative  showing  that  the  effect  of  bomidaries  is  to  retard  the  up- 

ward  migration  of  the  bubble.   For  x  <  -  ^  ,  i.e.  for  [i  >  2v    , 

d  -   becomes  positive,  while  d  t„  remains  negative.   Both   Id  , 
ool         ^       '         ol2  ^  '  ool 

and   |d  , ^1   tend  to  infinity  as  x  ->  -  1  .   Hence,  although  it  is 

OX  c_  '^\^ 

not  impossible  that  b  remains  negative  for  all  x  ,  it  will  cer- 
tainly tend  to  -co  for  x  ->  -1  ,  that  is  for  v  ->  0  ,  since 

Id  ,„|  ^>  o)  faster  than   Id  n  I  .   This  ca:i  be  seen  from  the  table 
'  ol2 '  '  ool' 

for  f  (x)   found  in  [6],  Fig.  1$,      We  can  conclude  therefore  that 
there  exists  a  critical  value  of  v   for  which  b<(2T)  '-•  0  ,  and  for 
V   smaller  than  this  critical  value  the  bubble  will  migrate  down- 
wards . 

3.     Shape  of  the  bubble. 

V7e  have  shovm  in  [2]  that  in  the  absence  of  boundaries  the 
bubble  would  flatten  in  the  first  approximation  and  would  become 
kidney  shaped  in  the  third  approximation.   This  is  because  "a^     <  0 
and  "a.  ^  >   0  .   The  effect  of  boundaries  on  the  shape  is  felt  only 
in  the  third  approximation  and  they  affect  only  the  value  of  a-,p, 
(see  equations  (7.2)).   By  inspecting  equations  (5.1^"^^)  it  nay  be 
verified  that  a-^p  >  0  for  sufficiently  small  v  .   Hence  ttp 


.;.=  •  '.i!- 


'3k- 


~  ^   ^22  ^    "^P  ^  * '  •  "^"^^  become  positive,  resulting  in  an  elongated 
bubble.   This  possibility  admits  of  a  very  simple  explanation.   The 
rigid  bottom  attracts  the  bubble  downvjards  while  a  free  surface  repels 
the  bubble,  also  downwards.   On  the  other  hand,  due  to  gravity,  the 
bubble  has  the  tendency  to  rise.   VJhen  the  bubble  is  near  the  bottom, 
the  upward  pull  due  to  gravity  is  larger  than  the  repelling  action 
of  the  free  surface.  As  a  result,  the  upper  part  of  the  bubble  will 
move  upward,  and  the  lower  part  -  dov/nward.   This  may  result  is  an 
elongation  of  the  bubble,  and  even  in  a  tear  into  tv/o  bubbles,  one 
below  the  other. 

Our  deductions  have  been  made  by  considering  only  the  .most 
significant  term.s  in  (7.2) .   Thus  in  discussing  periods,  we  con- 
sidered only  effects  due  to  first  order  terms,  and  neglected  to 
review  the  effect  of  second  order  and  third  order  tsrms.   In  study- 
ing the  migration,  we  also  neglected  the  third  order  terms.  All 
these  neglected  terms  may  play  a  significant  role.   However,  our 
conclusions  are  fully  confirraed  by  a  numerical  example  reported  on 
in  "Underwater  Explosion  Bubbles  IV." 


I 
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APPENDIX 


A,  Det©nnlnat:ion  of  the  ^, 


m 


(m+1) 


P  (cos  0).    i"      was  defined  in  Section 


Let   i  (r,e)  =  r'      . 

III  as  a  regular  harmonic  function  in  the  strip  -  u  <  ^  <  v  , 
satisfying  the  boundary  conditions 


(Al) 


h     +   i""  =  0  ,      for   ^  =  V   , 

'm    ^ra      *  ^      ' 


(h     +   b'"')^  =  0   ,        for   ^  =  -  u  , 

,  ^m   ^m '  ^     '  ^       ' 

V 


where. 


(A2) 


u 


^.  b, 


v  =  H  .  b 

CK 


^'"'  may  be  constructed  by  using  the  method  of  images.   Denote 
by  ^,   the  ^   coordinates  of  images  of  the  origin  in  the  planes 
^  =  V  ,  and  5  =  -  u  ,  lying  below  the  plane  z  =  -  u  ,  and  by 
"Z,    the   ^  coordinates  of  images  lying  above  the  plane   z  =  v  .   We 


have. 


(  -  2sw  ,   if  k  =  2s 


^ 


(A3) 


■\ 


<. 


V 


^k  - 


2sw  -  2u  ,    if  k  =  2s  +  1 

2sw  ,   if  k  =  2s 


2sw  +  2v  ,    if  k  =  2s  +  1 


3  =  0,  1  ... 
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where. 


(AU)        w  =  u  +  V  =  iitll   . 


It  is  easily  verified  that  (A3)  yields  the  t,     coordinates  of  all 
the  it^.ages  If  k  runs  from  1   on,  whereas  for  k  =  0,  ti,     =  T,     =0 

'  '  0     o 

Our  construction  Is  based  on  the  following  proposition:   If 
g(^,li,^)  =  g(^)   Is  a  regular  harmonic  function  everywhere  except  at 
the  origin,  then 


(A5)   G(5,h,^)  =  -  2g(^)  +  ^  (-l)''[g(^   ,-  ^)  +  g(^  -  ^   ) 

n=o 

+  SU  -   ?2n)  -  S(^2n+1-  ^^^ 

is  a  regular  harmonic  function  in  the  strip   -  u  <  ^  <  v  ,  and 
G(?,^,v)  =  -  g(v),  G,.(C,V^,u)  =  g^(-u)  .   We  specialize  g   to   (j)^(r,0) 
=  ^^{^,)^,K)    ,    obtaining 


OO  I 

Z 

n=o 


^m'^2n+l 


Here  the  ( * )   indicates  that  the  terms  involving  ^  =  ^  =  0  should 
L  ,  o    o 

be  omitted  from  the  series. 


Defining  r^^,  9^^,  r^^,  9j^  by 


.V^iU'il 


y    «> 
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(A7) 


i 


r^  +  kI  -  2r^^cos  0 
r^  +  Z^^  -  2r?j^cos  9 
(-l)^(r  cos  9  -  ^^) 


r,  cos  9,,  =  (-1)^{^  -  ^,J  =  (-Dd-  cos  9  -  ^,  )   , 


j^v^wo  ^^      -  ,-^,        ,^   -   ^J^ 


k' 


V 


we  now  get 


Substituting  In   (A6),  we  finally  obtain 


(A9) 


^:(r,9)  ^E*  (-l)^[r-(;t'^^.(-^  Wl)  -  -^^'^^.(-^  ^2n) 


n=o 


•^"an'^^Prn'^"-  «2n'  "  ^iirt"  V"=  «2n+l' ^  • 


B.    Expansion  of   .})"   near  r  =  0. 


m 


In  order  to  evaluate  (|)"(X,9)   as  a  power  series  in  o^  ,  an 


^m 


expansion  of  ^        in  powers  of  r  near  r  =  0  is  needed. 
By  definition  of  Legendre  Polynomials, 


CD 


^l-2:c  cos  9  +x^     ° 


x^P^(cos  9^)  ,     Ix|  <  1 


Hence,  on  replacing  x  by  —   , 

n 
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(Bl) 


/r.2.; 


r~-2xr^cos  6j^+x 


X  <  r 


n 


from  which  we  conclude  that 


MV®n^  -  k\^'d2^'^ 


vAA. 


r   -2xr  cos  9  +x 


x=0  . 


Since,  however,  by  (A7) 


r^  -  2xr  cos  0  +  x^  =  r^  -  2r(^^^  +  (-l)^x)cos  G 


we  get 


<p  (r  ,0  )  =  (-1)   — r(-sv-) 
^nr  n'  n    ^      rtil  ^  o^ 


n   V^ 


r  -2r^^cos  0  +K,^ 


Hence,  using  (Bl)  v;ith  r  replacing  x  ,  and  ^   replacing  r  , 
we  get,  on  carrying  out  the  differentiation. 


(E2'      M-n'«n'  =  l^nl"'  t   (-D'"*''"^'  ^^^   ^"'""''•^(-^  «' 


k=o 


-i  l«nl  • 


Similarly, 


(B3)       ^„(v«n)  =  l^nl-'  SJ-^'""-'"'"  ^ri^^;''^^'r^(-=  «'  . 


^  2  i?nl 


/.. 


-39- 


Substitution  In  (A9)  now  yields  the  desired  expansion. 


CO        1 
K=0 


where 


+  (_l)n+mi2^ _  (.i)n  j-^ji+l^,  ^  ^ 

1^2nl  l^2n+l' 


Using  the  definitions  of  ^  ,  ^  ,  we  get. 


> 
^"s^l 

CD 


/iDr-\  J        _    /      cT       xk+m+l    (k+m)l   U—    /    ^xS+m/,     ,     /    ,  v  k+in .    -  ( k+m+1) 


r:    (-l)^[(.l)^'+"^(s   +  ^)-(k+ra+l)_    (3   ^  £j^)-(k+rn+l)^ 


+ 
s=o 


^*  Expansion  of  d_^,   in  power  of  cf  . 

Since  v  +  ob  >  0  ,  \x  -    d]:>  >  0   ,    each  term  in  the  brace  of 
(B5)  may  bo  expanded  in  powers  of  <f  .      In  the  text  we  considered 
the  solution  for  X,    ^,    b  up  to  terras  in  cr    .   Since  the  expansion 
of  d .,   starts  with  a  term  in  cT  ^   ,  and  since,  in  view  of  the 
theorem  of  Section  III,  each  (j)^  is  multiplied  in  the  expression 
for   ij)  by   cf  to  at  least  m' th  power,  for  our  purposes  it  will 
be  sufficient  to  consider  the  expansion  of  the  brace  in  (B5)  up  to 
the  power  3  -  (2m+k+l)  .   Thus  we  only  need  to  consider  the  cases 
m=0,  k=0,  1,  or  2,  and  m  =  1,  k  =  0  . 


»>  •.  i  V 


-ko- 


Writing 


(CI) 


00 


mk 


:=o 


'^mk/  <^' 


,f 


we  have 


(C2) 


/ 


d   =  d  ,cr  +  d  -.cf^   +  . . . 
00    ool     oo3 

2 
^ol  =  ^ol?^  ■"  •  •  • 


V 


Uo   no 


^Ir,^'^        +  ... 


2  3 

Here,  the  terms  d  ^6       and  d  ,  ^(^       were  omitted,  since,  as  may  be 
'  oo2         ol3 

verified  they  are  zero.   We  now  get, 
f 


(C3) 


1 


d 


00 


-I  j       00       ,      ,   ^  S  CO 


-1 


-    (s   +  -#-)"l]    { 


00 


ol2 


•1  uu 

ij.(|J/TV)^    s=o  f^+^  1-L+v' 


CX) 


''         .,  1  iiL      /    n  ^s         CO 


-      (3     +-^)-^] 


^lo2  -   ^012 


\         ^003  =  2bid^3_2 


■  -*  >  /• 


-1^1- 


The  series  in  (C3)  can  be  expressed  in  terms  of  tabulated 
functions  of 


(CI;) 


^  V+|l  H+2q  * 


namely 


(C5)       f(x)  =21  (-1)    [(2S+1-X)   -"-(as+i+x)   -"] 
s=o 


=  2i;  [log    r   (ii2)  .  log   Ri?)]   -  f  tan  S| 


'dx 


^ 


and  its  derivatives.   We  find 


(C6) 


ool 


<   %12 


-  (f(x)  +  log  2)([J.  +  v) 

-  f'(x)(tx+v)"^ 

•  1^" 


-1 


d^23  =  -  ^T    ^^^  ■*"  fe  ^(3))(!^+v)'^ 


Here  ^  denotes  the  ^-Riemann  function. 

In  special  cases   (x  =  co ,  or  v  =  co  ,  formulae  (C6)  siuiplify 


to: 


r 


ool   2v 


(C7)     ^ 


V 


ol2 


'o23 


ool 


1 


8v- 


1_ 

21J. 


for  [i.  =  00  , 


(C8) 


<  d 


ol2 


'o23 


1 

8|? 


for  V  =  00  . 


-k2- 
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